
FINISH
/CLEARALL

*DO, it, 1, 5

/UIS, MSGPOP, 3		! Turns off note and warning messages
/UIS, ABORT, OFF	! Turns off status and cancellation messages
/KEYW, PR_SGUI, 1	! Prevent Solution done message - affects image plotting
/WINDOW, 1, OFF		! Turn graphics window off while running iterations
/PNUM, MAT, 1
/NUMBER, 1


max_elem = 4630		! Specify number of ceramic element in the model
n = max_elem + 2	! Starting index of the polarised elements
mesh_size = 21		! Model size
meshdivs = 21		! Model size
pole_volt = 20e5	! Apply voltage
W = 14853.44828		! Specify the coercive energy density

/INPUT, '2. start_script', 'txt', '', , 0
/INPUT, 'materials_SL_5', 'txt', '', , 0
/INPUT, '3. geom', 'txt', '', , 0
/INPUT, 'assign_SL_5', 'txt', '', , 0
/INPUT, '4. boundary_conditions', 'txt', '', , 0
/INPUT, '5. determine_res', txt, '', , 0
/INPUT, '6. write_res', txt, '', , 0
/INPUT, '7. shear_analysis', txt, '', , 0

/INPUT, '8. perm_analysis', txt, '', , 0

ESEL,U,MAT,,1  
/REPLOT 


FINISH
/CLEARALL

*ENDDO




